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FOREWORD 


This  report,  prepared  by  the  Lockheed-Callfornia  Company  under  Contract 
DOT-FA75WA-3707,  contains  a description  of  the  effort  in  which  program  KRASH 
was  used  to  model  and  analyze  a twin-engine,  low-wing  airplane  substructure 
subjected  to  a vertical  crash  Impact  condition.  The  test  was  performed  pre- 
viously by  NASA- Langley  as  part  of  a joint  FAA-NASA  effort  concerning  general 
aviation  airplane  crash  dynamics.  The  work  discussed  in  this  report  was 
administered  under  the  direction  of  the  Federal  Aviation  Administration  with 
H.  Spicer  acting  as  Technical  Monitor. 

Gil  Wittlin  of  the  Lockheed-Callfornia  Company  performed  the  analysis 
and  correlation  with  test  data.  Bill  LaBarge  of  the  Lockheed-Callfornia 
Company  assisted  In  the  initial  phases  of  the  math  model  development. 

Dr.  Robert  Hayduk  of  the  NASA-Langley  Impact  Dynamics  Research  Facility 
coordinated  the  effort  and  provided  structure,  test  and  film  data.  The  data 
presented  In  this  report  are  a partial  input  to  NASA  for  purposes  of  evaluat- 
ing the  state  of  the  art  in  light-airplane  crash  dynamics  modeling  and  analy- 
sis. The  Lockheed  effort  was  performed  under  the  supervision  of  J.E.  Wignot. 
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SUMMARY 


The  results  of  using  program  KRASH  to  analyze  the  dynamic  response  of  a 
twin-engine,  low-wing  airplane  substructure  subjected  to  a 27.5  ft/sec  verti- 
cal Impact  velocity  crash  condition  are  presented. 

Included  in  this  report  are  the  math  model  description,  pertinent  data 
from  the  test,  a comparison  of  analysis  versus  test  results  and  the  results 
of  a limited  parameter  sensitivity  study  using  program  KRASH. 

The  substructure  Is  modeled  symmetrically  in  KRASH  with  32  masses, 

57  member  elements  and  44  nonlinear  beam  element  degrees  of  freedom.  Floor 
and  occupant  pelvis  vertical  acceleration  responses  obtained  from  test  mea- 
surements are  compared  to  corresponding  analytical  results.  Occupant  chest, 
substructure  roof,  and  window-ledge  motions  are  also  compared.  Variations  in 
external  spring  (crushable  structure)  load-deflection  behavior,  seat  stiff- 
ness, occupant  representation,  analytical  filter  frequency  cutoff  and  model 
size  are  included  to  ascertain  their  effect  on  dynamic  behavior  for  the  par- 
ticular substructure  and  impact  condition  being  evaluated.  Conclusions  are 
presented  based  on  the  results  of  the  effort. 

Pertinent  math  model  and  computer  output  data  are  presented  in' Appendices 
A and  B. 
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SECTION  1 
INTRODUCTION 


Program  KRASH  has  been  validated  with  several  sets  of  full-scale  aircraft 
crash  test  data  (References  1 and  2).  The  availability  of  crash  test  data  for 
a twin-engine,  low-wing,  light-airplane  substructure  and  an  additional  impact 
condition  provided  another  opportunity  to  demonstrate  KRASH's  capability  to 
represent  the  significant  dynamic  response  characteristics  of  a survlvable 
crash  condition.  In  addition  to  KRASH,  the  twin-engine  airplane  substructure 
is  to  be  modeled  using  two  other  digital  computer  programs,  DYCAST  and  ACTION, 
designed  for  a structural  crash  dynamics  evaluation.  The  modeling  of  a par- 
ticular structure  for  a defined  Impact  condition  using  three  different  cur- 
rent crash  dynamics  computer  programs  provides  an  opportunity  to  compare  the 
requirements  for:  model  size,  input  data,  ease  of  modeling,  output  data, 
analysis  versus  test  results,  machine  costs  and  machine  time.  This  report 
concerns  itself  solely  with  the  KRASH  modeling  results  and  comparison  with 
the  substructure  test  results. 
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SECTION 2 

MODEL  AND  TEST  DATA 


2.1  KRASH  MODEL 

The  subsCructure  representing  F.S.  135  to  F.S.  181  for  a typical 
twin-engine  low-wing  airplane  (Figure  2-1)  is  shown  in  Figure  2-2. 

The  substructure  including  occupants  is  modeled  symmetrically  (about 
B.L.  0.0)  in  program  KRASH.  The  KRASH  substructure  and  occupant  representa- 
tions are  shown  in  Figures  2-3  and  2-4.  The  floor  masses  at  locations  2 
through  6 and  8 through  12  have  external  springs  to  represent  ground  contact 
and  crushing  of  the  underside  structure.  Only  half  the  structure  is  shown  in 
Figure  2-3  since  the  model  and  Impact  conditions  are  treated  symmetrically. 

The  dash  lines  jutting  laterally  from  14  and  18  are  tension-only  members  to 
represent  the  lateral  tie  rods  (see  Figures  2-5  and  2-6).  For  clarity, 
tension-only  members  representing  occupant  seat  belts  connecting  mass  5 to 
mass  30  and  mass  11  to  mass  30  are  not  shown  in  Figure  2-4.  Compression- 
only  members  representing  the  seat  cushion  and  pan  connect  mass  30  to  mass 
29.  Masses  29,  30,  31  and  32  represent  the  seat,  occupant  lower  torso, 
occupant  upper  torso,  and  DRI,  respectively.  Details  of  the  modeling  including 
sample  calculations  are  shown  in  Appendix  A.  Computer  program  input  and  output 
data  are  provided  in  Appendix  B. 

• KRASH  Substructure  model  size: 

Total  number  of  masses  - 32 

Total  number  of  beam  elements  - 57 

Number  of  nonlinear  degrees  of  - 44 

freedom  associated  with  beam 
elements 
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FRONT  SPAR  ATTACH  FACE  RE  AN  SPAA  ATTACH  FACE 

MAIN  SPAN  DATUM  LINE 


SUBSTRUCTURE 

REGION 


Figure  2-1.  - Twin-engine,  low-wing  airplane. 


Figure  2-2.  - Lower  fuselage  structure  for  twin-engine 
low-wing  airplane  (F.S.  135  - F. S.  181). 
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Figure  2-3.  - KRASH  fuselage  model 


Figure  2-4.  - KRASH  floor-seat-occupant  model 
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Velocity  ■ 27.5  ft/sec  vertical 

Pitch  angle  ■ 3/4°  pitch-up 

Yaw-angle  " 0° 

Roll  angle  - 0° 

i Table  2-1  shows  the  location  of  the  test  and  analysis  data  points. 

j • 

The  nonlinear  deflection  values  used  as  input  into  the  KRASH  model 
were  obtained  directly  from  KRASH  calculations.  The  deflection  values  are 
contained  as  part  of  the  'Model  Parameter  Data'  print,  provided  in  Appendix  B. 

2.2  TEST  DATA 

The  substructure  crash  test  data  and  film  analyses  were  provided  by  the 
NASA-Langley  Impact  Dynamics  Research  Facility.  Figures  2-5  and  2-6  show  the 
front  and  rear  views,  respectively,  of  the  post-crash  test  condition  of  the 
substructure.  Figures  2-7  and  2-8  show  the  floor  structure  and  a portion  of 
the  roof  structure,  respectively.  The  summary  of  NASA-provided  test  data  is 
presented  in  Table  2-2.  Figure  2-9  shows  the  motion  histories  obtained  from 
high-speed  film  analysis,  performed  by  NASA-Langley,  for  the  occupant  chest, 
roof,  and  window-ledge  locations.  Figure  2-10  shows  the  reduced  test  data 
accelerometer  histories  for  four  floor  locations  and  the  two  occupant  (left 
and  right  side)  pelvis  locations.  The  channels  record  D.C.  data.  The  reduced 
data  are  obtained  by  NASA  using  a least-square  fit  (LSF)  filtering  technique. 
All  the  test  data  are  for  vertical  response  channels  and  motions  since  this 
direction  represents  the  predominant  mode  of  response  for  this  type  of  impact 
condition.  Correspondingly,  all  correlation  with  analysis  is  limited  to 
responses  in  the  vertical  direction. 

, The  substructure  test  was  previously  performed  by  NASA  as  part  of  a Joint 

general  aviation  crash  test  program  initiated  by  the  FAA  and  NASA.  The  impact 
dynamics  test  facility  and  test  procedures  are  described  in  Reference  3. 


LOCKHEED 

m.inlNit  f n M P A N V 


2-5 


TABLE  2-1.  - LOCATION  OF  TEST  AND  ANALYSIS  DATA  POINTS 


Teat  Data  Points 


Channel 

Location 

mm 

W.L. 

B.L. 

A1 

Floor 

140. 

-16. 

6L 

A2 

Floor 

151. 

-16. 

24L 

B2 

Floor 

151. 

-16. 

6L 

D1 

Floor 

162. 

-16. 

24L 

D9 

Left  Pelvic 

- 

- 

- 

DIO 

Right  Pelvic 

- 

- 

- 

(a) 

Analysis  Data  Points 

Mass 

Location 

F.S. 

W.L. 

B.L. 

2 

Floor 

140. 

-16. 

6L 

, 

3 

Floor 

151. 

-16. 

6L 

9 

Floor 

151. 

-16. 

20L 

10 

Floor 

163. 

-16. 

20L 

30 

Lower  Torso 

- 

- 

- 

^Analysis  is  performed  with  symmetrical  half  model,  thus  left 
side  ■ right  side. 
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2-5.  - Post-crash  condition  of  Figure  2-6.  - Post-crash  condition  of 

substructure,  front  view  substructure,  rear  view 


Figure  2-7.  - Post-crash  condition  of  substructure,  floor 


Figure  2-8.  - Post-crash  condition  of  substructure,  partial  of  roof 


TABLE  2-2.  - SUMMARY  OF  SUBSTRUCTURE  TEST  DATA 


Magnitude 

Impact  velocity^ 

330  ln/sec;  27.5  ft/sec 

Rebound  height  ^ 

5 inches 

Rebound  velocity ^ 

62.1  ln/sec;  5.2  ft/sec 

Duration  of  ground  con  tact ^ 

0.058  + 0.0015  sec. 

Attitude  at  impact 

0.75  degrees  pitch-up 

Window  ledge:  peak  motion 

1.64  in. 

final  deflection^6) 

0.65  in. 

(a) 

Roof  center:  peak  motion 

3.65  in. 

final  deflection^6) 

0.99  in. 

Plots <a) 

- Occupant  chest  motion  ve.  time 

- Roof  motion  vs.  time 

- Window-ledge  motion  vs.  time 

Accelerometer  traces  (kMO 

Floor  structure  - Al,  A2,  B2,  Dl 

Occupant  pelvis  - D9,  DIO 

Obtained  from  high  speed  film  analysis 

^k^Least  Square  Fit  (LSF)  filtered  data 
( c ) 

D.C.  Accelerometers 

^^No  significant  roll  or  yaw  motion  obtained  from  the  film  analysis 

(e) 

Post-test  measurement 
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TIME,  SEC 

Figure  2-9.  - Motion  histories  obtained  from 
high-speed  film  analysis 
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ACCELERATION.  G ACCELERATION.  G ACCELERATION.  G ACCELERATION.  G 


(1)  -84.9  9 CHANNEL  A1 


130  b-  (2) 


(2)  69.0  g (4)  40.9  g CHANNEL A2 


i i i i i i i i i 


I I I I t I 


0 0.02  a04  0.06  0.08  0. 


120  fa  «« 


TIME,  SECONDS 


(1)  -116.7  g (3)  26.9  g (6)  67.6  g 


(2)  17-4 g (4)  —81.0 g CHANNEL  B2 

I I I I I I I I I I I I I I I I 


0 0.02  0.04  a 06  0.08  0. 


TIME.  SECONDS 


(2)  (4) 


(1)  -80.8  g (3)  -24.2  g (6)  -10.2  g 


(2)  28.7  g (41  17.4  g CHANNEL  01 


TIME.  SECONDS 


(1)  -464 g (3)  -44.7 g CHANNEL  D9 


<2>  -31.1  g 

(1)  —62.7 g (3)  -43.6g  CHANNEL  D10 


(21  -20.4  g 


0.02  0.04  aoe  aos  o.io  0.12  o.u  o.ie  0.18  0.20  0.22 

TIME,  SECONDS 


Figure  2-10.  - Accelerometer  traces 
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SECTION  3 

KRASH  RESULTS  AND  COMPARISON  WITH  TEST  DATA 


3.1  ACCELERATIONS 

Table  3-1  shows  a comparison  of  the  analysis  and  test  peak  accelerations 
for  the  occupant  pelvis  and  floor  responses.  Figures  3-1,  3-2  and  3-3  pre- 
sent the  analysis  and  test  response  histories  for  the  occupant  pelvis  and 
floor  locations.  The  data  in  Table  3-1  show  that  the  analytically  obtained 
occupant  pelvis  peak  accelerations  agree  within  19  percent  for  the  initial 
peak  and  6.8  percent  for  the  second  peak,  when  compared  to  the  average 
of  the  left  and  right  side  occupant  test  measured  responses.  It  is  noteworthy 
that  the  left  and  right  side  measured  test  responses  for  a "symmetrical" 
impact  condition  differ  by  16  percent  and  2.5  percent,  respectively,  for  the 
primary  and  secondary  peaks.  The  difference  between  test  results  is  nearly 
the  same  magnitude  as  the  difference  between  the  analysis  and  test  results. 

The  analytical  data  (Figure  3-1)  show  the  same  "camel  hump"  response  phen- 
omenon as  that  which  is  exhibited  in  the  test  data.  In  addition  the  two  peaks 
obtained  by  analysis  occur  within  4 to  6 milliseconds  of  the  time  at  which 
they  occur  during  the  test.  The  test  and  analysis  results  show  a decay  to  zero 
acceleration  after  the  secondary  peak  is  reached.  The  decay  rate  is  approxi- 
mately the  same  for  the  analysis  and  test  data,  with  the  analytical  data  pre- 
ceding the  test  data  by  several  milliseconds.  The  analytical  results  for  the 
occupant  pelvis  response  presented  in  Table  3-1  and  Figure  3-1  are  based  on  a 
low-pass  filter  frequency  cutoff  of  100  Hz  as  noted  in  Section  2.  The  use  of 
a higher  frequency  cutoff  value  will  result  in  analytical  results  which  will 
be  higher  and  consequently  closer  to  the  test  results.  For  example,  a 150  Hz 
frequency  cutoff  in  the  analysis  will  increase  the  primary  and  secondary  occu- 
pant response  by  5 percent  and  2 percent,  respectively.  This  change  will 
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TABLE  3-1.  COMPARISON  OF  ANALYSIS  AND  TEST  PEAK  ACCELERATIONS 


Pelvis 


1st  Peak 
2nd  Peak 


Pelvis 


1st  Peak 
2nd  Peak 

Pelvis 

1st  Peak 
2nd  Peak 

Floor 

Forward  Inboard 


Forward  Outboard 


Forward  Inboard 


Aft  Outboard 


Test 

Accel  (Time)^8) 

Analysis  , . 
Accel  (Timer  ; 

Percent (b) 
Differences  (Time) 

(D9)(d) 

(Mass  30)(e) 

-45.4  (.018) 

-39.7  (.017) 

-12.6  (.001) 

-44.7  (.030) 

-41.1  (.026) 

- 8.1  (.004) 

(D10) 

(Mass  30) 

-52.7  (.020) 

-39.7  (.017) 

-24.7  (.003) 

-43.6  (.032) 

-41.1  (.026) 

- 5.7  (.006) 

(Avg.  of  D9  6 D10) 

(Mass  30) 

-49.0  (.018-. 020) 

-39.7  (.017) 

-19.0  (.001-. 003) 

-44.1  (.030-. 032) 

-41.1  (.026) 

- 6.8  (.004-. 006) 

(Al) 

(Mass  2) 

-84.9  (.004) 

-89.0  (.004) 

4.8  (0) 

9.3  (.010) 

43.9  (.028) 

(O 

(A2) 

(Mass  9) 

-128.7  (.010) 

-101.0  (.011) 

-21.5  (.001) 

69.0  (.018) 

61.8  (.025) 

-10.4  (.007) 

(B2) 

(Mass  3) 

-115.7  (.004) 

-90.7  (.004) 

-21.6  (0) 

25.9  (.032) 

34.9  (.026) 

34.0  (.006) 

(Dl) 

(Mass  10) 

-80.8  (.008) 

-99.2  (.007) 

22.8  (.001) 

28.7  (.017) 

37.1  (.018)<f) 

29.3  (.002) 

3 ) 

Accelerations  In  G's,  time  in  seconds  after  impact 

b)Percent  difference  - (ana ly8is  v«lue:te8t  value)X 

\ test  value  ' 

c ) 

Test  trace  shows  no  response  after  .012  seconds 


'Test  channel  numbers  in  parenthesis 
^Analysis  mass  numbers  in  parenthesis 
^Peak  at  .026  msec  * 53.7  g 

NOTE:  Difference  between  test  results  for  right  (DIO)  versus  left  (D9) 
side  *16.1%  (1st  peak)  and  2.5%  (2nd  peak). 


3-2 


ACCELERATION  G 


0.02  0.04  0.06 

TIME,  SECONDS 


Figure  3-1.  - Comparison  of  occupant  pelvia  vertical 
acceleration,  test  versus  analysis. 

result  in  closer  correlation  agreement  with  test  data.  A discussion  of  the 
sensitivity  of  analysis  results  to  the  selection  of  an  analytical  filter  cut- 
off frequency  is  discussed  in  Section  4.4. 

From  Table  3-1  it  can  be  seen  that  the  analysis  initial  peak  acceleration 
values  agree  within  approximately  ±22  percent  with  measured  test  data  obtained 
from  the  four  floor  dc  channel  accelerometer  locations  which  were  chosen  for 
comparison  by  NASA- Langley . NASA  considers  dc  accelerometers  to  be  overall 
more  reliable  than  ac  accelerometers,  particularly  the  behavior  after  the 
initial  peak  response.  The  time  of  occurrence  of  the  initial  peak  values 
obtained  by  analysis  agrees  with  the  test  data  peak  value  occurrences  within 
one  (1)  millisecond.  Figures  3-2  and  3-3  show  the  similarity  in  response 
shapes  for  the  test  and  analysis  results,  particularly  in  the  initial  response 
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Figure  3-2.  - Comparison  of  fuselage  floor  vertical 
acceleration,  test  versus  analysis. 
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Figure  3-3.  - Comparison  of  fuselage  floor  vertical 
acceleration,  test  versus  analysis 
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regime.  The  data  presented  in  Table  3-2  show  that  the  secondary  (rebound) 
peak  accelerations  and  times  of  occurrence  of  these  peak  values,  obtained  by 
analysis,  agree  within  -10  to  +34  percent  and  within  2 to  7 milliseconds, 
respectively,  with  the  corresponding  measured  test  data.  The  analytical  floor 
results  provided  in  Table  3-1  and  Figures  3-2  and  3-3,  like  the  occupant  pelvis 
data,  are  based  on  a 100-Hz,  low-pass  filter  cutoff  frequency.  For  a 150-Hz 
cutoff  frequency,  the  analytical  results  will  be  higher.  However,  the  compari- 
sons for  the  Initial  peak  values  show  the  analysis  results  to  be  within  -11.7 
to  +34  percent  of  the  corresponding  test  values.  While  the  test  ac  acceler- 
ometers are  not  used  in  the  base  case  comparisons  presented  in  Table  3-1  and 
Figures  3-1,  3-2,  and  3-3,  it  is  of  significance  to  note  that  the  measured 
test  values  vary  by  as  much  as  16  to  55  percent  between  the  left  and  right 
sides  of  the  substructure  when  both  ac  and  dc  accelerometers  are  included  in 
the  data.  This  difference  is  equal  or  greater  than  the  difference  between 
analysis  and  test  results.  The  effect  of  cutoff  frequency  on  floor  responses 
is  discussed  in  Section  4.4. 

3.2  DISPLACEMENTS 

Table  3-2  and  Figure  3-4  show  the  comparison  of  analysis  and  test  motions 
for  the  occupant  chest,  the  substructure  roof,  and  the  substructure  window- 
ledge.  The  analytical  results  are  based  on  obtaining  the  vertical  displace- 
ment of  mass  31  (Figure  2-4)  for  the  chest  motion,  the  average  of  the  vertical 
displacement  of  masses  24  through  28  (Figure  2-3)  for  the  substructure  roof 
motion,  and  the  average  of  the  vertical  displacements  of  masses  14  through  17 
(Figure  2-3)  for  the  substructure  window-ledge  displacements.  The  peak  values 
obtained  by  analysis  are  within  -13  percent,  -27.4  percent  and  +43.3  percent 
of  the  peak  displacements  for  the  chest,  roof,  and  window-ledge  motions, 
respectively,  obtained  from  analysis  of  high-speed  films  of  the  crash  test. 

The  time  of  occurrence  of  the  peak  values  obtained  by  analysis  agree  with  the 
film- analyzed  test  peak  occurrences  within  1 to  3 milliseconds.  The  motion 
histories  (Figure  3-4)  depict  overall  agreement  between  analysis  determined 
and  test-observed  trends,  even  though  the  peak  values  show  disagreements  of 
between  13  and  43  percent. 
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TABLE  3-2.  - COMPARISON  OF  ANALYSIS  AND  TEST  PEAK  MOTIONS 


Test 

Displacement 

(time) 

Analysis 

Displacement 

(time) 

Percent 

Difference 

(time) 

Occupant  Chest 

8.40  (.030) 

7.30  (.033) 

-13.1  (.003) 

Roof 

3.65  (.015) 

2.65  (.013) 

-27.4  (.002) 

Window  Ledge 

1.64  (.01-. 015) 

2.35  (.012) 

43.3  (.002-. 003) 

(a) 

Displacement  in  inches,  time  in  seconds  after  impact 

(b) „  _ ,,,,  analysis  value  - test  value  ,,  ,,,- 

Percent  difference  - i — - — = X 100 

test  value 


The  analytical  results  indicate  that  lower  fuselage  deformation  of  up  to 
1.4  inches  occurs,  which  appears  to  be  approximately  twice  as  high  as  the 
average  deformation  observed  in  the  test  high-speed  film  review.  A review  of 
the  post-crash  test  condition  of  the  structure  Indicates  that  overall  defor- 
mation of  the  structure  may  not  have  exceeded  0.6  inches  based  on  the  average 
of  the  front  and  rear  deformation. 

The  analysis  shows  that  all  the  elements  of  the  structure  behave  lin- 
early, except  for  the  seat  cushion  and  pan,  the  fuselage  underside  and  several 
elements  of  the  shell  structure  (beams  7-13,  8-14,  9-15,  10-16,  11-17  and 
12-18,  Figure  2-3).  The  postcrash  condition  of  the  substructure  indicates 
little  in  the  way  of  permanent  structural  deformation.  The  analysis  results 
show  a Dynamic  Response  Index  (DRI)  value  of  28.  A DRI  magnitude  of  28  indi- 
cates a high  potential  (»50  percent)  for  a spinal  vertebrae  comp  res  si  on- type 
injury  to  occur.  The  available  test  data  do  not  allow  for  an  assessment  of 
this  type  of  injury  potential. 

Variation  of  the  analytical  filter  cutoff  frequency  has  no  effect  on  the 
motion  results  since  filtering  is  a post-processing  procedure  in  KRASH  and 
does  not  alter  any  of  the  basic  results. 
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3.3  COMPUTER  COST  AND  PERFORMANCES 

The  KRASH  analysis  was  performed  on  an  IBM  370-3330  digital  conputer. 

For  the  size  case  and  impact  condition  analyzed  (see  Section  2.0)  the  approxi- 
mate cos t per  computer  run  is  $135.  The  analysis  (machine)  time,  including 
all  printed  and  plotted  data,  is  6.9  minutes  per  run. 
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KRASH  SENSITIVITY  STUDY  RESULTS 

An  analysis  using  program  KRASH  Is  based  on  the  premises  that: 

• Nonlinear  load-deflection  behavior  can  be  approximated 

• Only  a portion  of  the  major  structural  elements  need  be  modeled 
nonlinearly  for  post-failure  behavior. 

Since  the  use  of  KRASH  Is  designed  to  obtain  overall  responses,  depict 
significant  dynamic  phenomenon,  and  represent  response  trends  using  approxi- 
mate modeling  techniques,  it  is  understandable  that  different  users  might 
model  aircraft  structure  or  select  input  data  which  varies  in  some  respects. 

To  determine  the  consequence,  insofar  as  dynamic  response  results  are  con- 
cerned, of  establishing  a model  and  an  analysis  which  might  vary,  a limited 
sensitivity  investigation  was  performed.  The  base  case  results  which  best 
represent  the  estimated  properties  of  the  structure  and  occupant  are  described 
in  Section  3.0.  In  this  section  the  results  of  varying  the  following  param- 
eters are  discussed: 

• External  spring  load-deflection  representation 

• Occupant  axial  stiffness 

• Seat  maximum  force  level 

• Analytical  filter  cutoff  frequency  selection 

• Math  model  size 

There  was  little  in  the  way  of  nonlinear  behavior,  other  than  in  the 
crushable  lower  fuselage,  occupant  seat  cushion  and  pan  deflection,  and  several 
elements  in  the  airframe  shell,  noted  in  the  analysis  of  the  substructure  for 
the  defined  impact  condition.  Consequently,  the  above  noted  parameters  are 
considered  representative  of  areas  wherein  different  users  may  have  differences 
of  opinion  with  regard  to  modeling  representations. 
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4.1  EXTERNAL  SPRING  LOAD-DEFLECTION  REPRESENTATION 


The  representation  of  the  crushable  lower  fuselage  structure  within  the 
framework  of  KRASH's  external  spring  Input  requirements  Is  perhaps  the  most 
sensitive  and  critical  concern  for  modeling.  The  data  for  the  base  case 
representation  were  obtained  following  the  procedures  outlined  in  References  4 
and  5.  A sample  calculation  is  shown  in  Appendix  A.  The  most  likely  area  of 
variation  in  a user's  thinking  would  be  in  the  representation  of  the  slope  of 
the  load-deflection  curve  in  the  nonlinear  region.  The  base  case  assumes  a 
constant  load  in  the  nonlinear  region  until  bottoming  occurs.  This  curve 
(condition  A)  is  shown  in  Figure  4-1,  along  with  possible  user  input  varia- 
tions. Condition  B provides  for  a constant  20-percent  increase  in  the  assumed 
maximum  forces  associated  with  the  crushable  structure.  Conditions  C and  D 
represent  sloping  20-percent  force  changes  (increasing  and  decreasing,  respec- 
tively) up  until  bottoming  occurs.  All  the  curves  are  based  on  a 0.1-inch 
initial  linear  region  and  a 3-inch  deflection  before  bottoming  occurs.  Both 
values  are  determined  following  the  procedure  outlined  in  References  4 and  5. 
The  results  would  not  be  altered  to  any  significant  degree  if  the  linear 
deflection  value  were  in  the  range  of  0.01  to  as  much  as  0.3  inch,  which 
leaves  significant  margin  for  modeling  variation.  Since  the  maximum  deflec- 
tion of  the  crushable  structure  in  the  base  analysis  does  not  exceed  1.6  Inches 
the  selection  of  3 inches  for  bottoming  to  occur,  is  also  of  little  conse- 
quence in  this  analysis.  Thus,  it  is  reasonable  to  assume  that  the  selection 
of  the  peak  force  and  the  rate  at  which  force  varies  with  deflection  will  be 
most  influential  on  the  results  for  this  analysis.  The  results  of  this  param- 
eter sensitivity  investigation  are  shown  in  Figure  4-2.  The  20-percent  con- 
stant increase  in  external  spring  force  maximum  load  results  in  an  approxi- 
mately 15  to  20-percent  increase  in  peak  floor  accelerations  at  the  four  floor 
dc  channel  accelerometer  locations.  Changing  the  slope  (conditions  C and  D) 
has  a substantially  less  effect  on  the  peak  responses.  The  sloping  ±20  per- 
cent change  results  in  a variation  of  the  peak  floor  accelerations  of  be- 
tween ±3  to  ±11  percent.  The  occupant  pelvis  responses  are  not  affected 
by  these  changes  because  in  this  model  and  under  the  defined  impact  conditions 
the  pelvis  responses  are  Influenced  primarily  by  the  seat  and  occupant 
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Figure  4-1.  - Normalized  force  versus  external  spring 
load-deflection  characteristics. 

characteristics.  As  can  be  anticipated,  the  maximum  deflection  of  the  external 
spring  decreases  as  the  force  level  increases  and  vice  versa.  For  the  range 
of  external  spring  characteristics  investigated,  the  average  floor  deformation 
ranges  from  1 to  1.5  inches.  The  range  of  test-measured  peak  accelerations  is 
also  shown  in  Figure  4-2.  The  analysis  results  for  the  range  of  load-deflection 
behavior  investigated  are  within  the  range  of  the  measured  ac  and  dc  acceler- 
ometer data.  When  considering  dc  accelerometer  data  alone  the  analysis  peak 
values  are  no  more  than  30  percent  different  than  the  test  peak  values  except 
at  the  aft  outboard  floor  location  where  the  difference  is  50  percent  and  the 
ac  and  dc  accelerometers  differ  by  75  percent  (using  dc  accelerometer  values 
as  base  values) . The  results  presented  in  Figure  4-2  Indicate  that  the  sensi- 
tivity of  the  analysis  results  are  comparable  to  the  sensitivity  of  the  test 
measurements.  The  results  further  indicate  that  all  four  external  spring 
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The  representation  of  the  occupant  ia  controlled  by  the  user  In  program 
KRASH  by  the  selection  of  mass,  area,  and  damping  properties  assigned  to  the 
corresponding  masses  and  beam  elements.  The  lower  and  upper  torso  masses 
(30  to  31  in  Figure  2-4)  are  assigned  mass  weight  and  inertia  properties  as 
described  in  References  5 and  6 and  in  Appendix  A.  The  connecting  link  between 
masses  31  and  32  represents  the  properties  of  the  upper  spine.  The  axial 
stiffness  of  this  element  is  determined  from  the  area,  modulus  of  elasticity, 
and  length  properties  assigned  to  this  element.  The  data  in  Reference  7 
suggest  that  the  human  upper  torso  has  specific  frequency  and  damping  values. 


LOCKHEED 

C»U»0"NI»  COMPANY 


Trans  laced  Into  KRASH  data  this  requires  the  axial  frequency  and  the  damping  of 
the  human  spine  to  be  14  Hz  and  41  percent  of  critical,  respectively.  Depend- 
ing on  the  actual  properties  of  the  anthropomorphic  dummy  used  In  a particular 
test  It  Is  possible  that  the  frequency  and/or  the  damping  can  be  different 
than  the  values  suggested  for  humans.  In  this  particular  case  the  spinal 
axial  stiffness  is  not  defined  for  the  dummies  used  In  the  test. 

For  the  human  properties  defined  In  Reference  7,  KRASH  Is  coded  to  com- 
pute the  desired  frequency  and  damping.  For  all  other  spinal  stiffness,  the 
user  can  easily  determine  the  proper  values  to  input.  Knowing  the  desired 
frequency,  the  input  mass  values,  and  the  length  of  the  spinal  element, 

combinations  of  modulus  of  elasticity  and  cross  section  area  may  be  determined 
for  input  into  KRASH.  Damping,  as  a percent  of  critical,  can  be  user  selected 

for  any  individual  element. 

To  evaluate  the  consequence  of  different  spinal  axial  properties  on  the 
response  of  the  structure  and  occupant,  a variation  in  upper  torso  axial  fre- 
quency from  14  to  61  Hz  was  examined.  The  results  of  this  investigation  are 
presented  in  Figure  4-3.  As  shown  in  Figure  4-3,  as  the  axial  frequency 
increases  the  analytical  results  tend  to  reproduce  the  two  peaks  (camel 
hump)  effect  exhibited  in  the  test  data  for  the  occupant  lower  torso.  Using 
an  axial  frequency,  of  61  Hz.,  a maximum  seat  force  of  5570  pounds  and  0.41 
damping  the  analysis  results  are  20  to  30  percent  lower  than  the  test  data. 

For  a higher  seat  force,  the  occupant  lower  torso  acceleration  values  increase. 
For  a lower  axial  stiffness  the  occupant  lower  torso  initial  peak  value 
increases  while  the  second  peak  value  decreases,  until  it  actually  goes  negative 
at  14  Hz.  Changing  the  damping  value  while  holding  all  other  parameters  con- 
stant tends  to  affect  the  second  peak  value  more  than  the  initial  peak  value. 
During  a separate  study  using  a simplified  5-mass  system,  changing  damping 
from  0.31  to  0.04  increased  the  occupant  lower  torso  first  peak  value  by  2.5 
percent  and  the  second  peak  value  by  29  percent. 

The  parameter  sensitivity  analyses  involving  the  occupant  representation 
illustrates  the  importance  of  providing  input  data  representative  of  the  ele- 
ments being  modeled.  It  may  well  be  that  the  properties  of  test  dummies 
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Figure  4-3.  - Lower  torso  (pelvis)  response  as  a function 
of  occupant  axial  stiffness 

differ  from  the  corresponding  human  properties.  It  appears  from  the  results  of 
this  investigation  that  the  "camel  hump"  effect  noted  in  the  test  data  means 
that  the  dummy  spine  is  stiffer  than  originally  anticipated  in  setting  up  a 
representation  based  on  Reference  7 data.  The  fact  that  the  occupant  upper 
torso  responses  follow  very  closely  the  occupant  lower  torso  responses  in  both 
magnitude  and  time  of  occurrence,  supports  the  contention  that  the  spinal 
connection  of  the  dummies  used  in  the  substructure  test  are  relatively  stiff. 

4.3  SEAT  MAXIMUM  FORCE 

The  seat-cushion  and  pan-stiffness  properties  used  in  the  analysis  are 
representative  of  those  installed  in  the  particular  type  of  airplane  used  in 
the  substructure  test.  Figure  4-4  shows  a general  force-deflection  curve 
for  a light  aircraft  passenger  seat.  Since  the  curve  provided  in  Figure  4-4 
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Figure  4-4.  - Typical  force-deflection  curve  for 
a light-aircraft  passenger  seat. 


is  not  the  actual  load-deflection  for  the  seat  installation  in  the  substructure 
test,  an  investigation  into  the  effect  of  the  variation  in  maximum  seat  force  on 
occupant  response  was  performed.  The  change  in  cushion  and  pan  stiffness  ef- 
fects the  occupant  response  to  a lesser  extent  than  the  force  cutoff  value  in 
this  particular  situation  and  was  not  fully  explored.  Figure  4-5  shows  the 
results  of  the  force  cutoff  changes.  The  occupant  pelvis  responses  increase 
as  the  maximum  force  level  increases.  For  a maximum  cutoff  force  of  5570 
pounds  the  analysis  results  are  approximately  30  to  19  percent  lower  than  the 
average  of  the  measured  responses  for  the  two  peak  accelerations.  For  the 
6455  pounds  used  in  the  base  analysis  the  analytical  results  are  approximately 
19  and  7 percent  lower  than  the  average  of  the  measured  responses  for  the  two 
peak  accelerations.  If  the  cutoff  force  was  as  high  as  9000  pounds,  the 
analytical  results  would  be  8 and  18  percent  higher  than  the  average  of  the 
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OCCUPANT  DAMPING  - 41  PERCENT  OF  CRITICAL 
OCCUPANT  AXIAL  FREQUENCY  - 61  HZ 


ANALYSIS  FILTER  ft  100  HZ  LOW-PASS 
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Figure  4-5.  - Lower  torso  (pelvis)  response  as  a function  of 
seat  maximum  force  level. 

measured  responses.  The  ranges  investigated  represent  a deviation  of  approxi- 
mately -15  to  +28  percent  from  the  nominal  6455  pounds  used.  For  the  three 
force  levels  evaluated  (557C,  6455,  and  9000  pounds),  the  initial  peak  accel- 
eration level  obtained  by  analysis  for  the  occupant  upper  torso  is  approxi- 
mately 6 percent,  17.4  percent,  and  58  percent  higher,  respectively,  than  the 
corresponding  average  of  the  head  and  chest  measured  values.  The  upper  torso 
response  is  more  akin  to  the  DRI  value  than  the  lower  torso  response.  Con- 
sequently, the  analysis  may  show  a more  conservative  DRI  value  than  is 
indicated  by  test  response  data.  User-selected  values  of  between  5500  and 
6500  pounds  appear  to  provide  reasonable  representations  of  the  seat  load- 
deflection  characteristics. 

4.4  ANALYTICAL  FILTER  CUTOFF  FREQUENCY 

Analysis  of  the  current  substructure  test  data  indicated  that  a 100-Hz, 
low-pass  digital  filter  is  approximately  equivalent  to  the  least-square-fit 
(LSF)  filtering  performed  by  NASA-Langley . The  LSF  data  were  used  for  the  tabu- 
lation of  test  results.  Previous  NASA  evaluation  of  equivalent  filtering 
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indicated  that  180  Ha,  low-pass  filtering  could  be  equivalent  to  the  LSF  data. 
Recognizing  that  equivalent  filtering  is  difficult  to  define  because  of  the 
different  data  reduction  processes  that  are  used  and  the  wide  range  of  types 
of  structures  or  elements  being  considered,  an  evaluation  was  performed  to 
ascertain  the  effects  of  the  analytical  results  and  subsequent  correlation  with 
test  results  of  different  analytical  filter  cutoff  frequencies.  As  noted 
earlier,  KRASH  filtering  is  a post-processing  technique  and  doesn't  change  any 
computed  values. 

Figure  4-6  shows  the  result  of  the  analysis  for  a cutoff  frequency  range 
from  100  to  150  Hz.  For  the  floor  responses  the  analytically  obtained  peak 
values  increase  by  approximately  4 to  18  percent  in  changing  from  a 100-Hz  to 
150-Hz  cutoff  frequency.  However,  as  can  be  observed  in  Figure  4-6,  the 
results  are  for  the  most  part  still  within  the  range  of  recorded  test  values 
at  the  respective  locations.  The  analytical  responses  for  the  occupant  lower 
torso  are  less  sensitive  to  the  cutoff  frequency  change.  The  data  in  Figure 
4-6  shows  a variation  of  less  than  5 percent  for  the  initial  lower  torso  peak 
value  and  less  than  2 percent  for  the  second  peak  value.  The  occupant  re- 
sponses exhibit  lower  frequency  (broader  response)  characteristics  than  the 
floor  responses  and  consequently  are  not  expected  to  be  as  sensitive  to  the 
higher  cutoff  frequencies. 


4.5  MODEL  SIZE  VARIATION 

The  base  32  mass,  57  member  symmetrical  math  model  shown  in  Figures  2-3 
and  2-4  provides  an  adequate  representation  of  the  substructure  and  Impact 
condition  evaluated,  as  attested  to  by  the  close  agreement  with  available  test 
data.  The  results,  analytical  and  test,  Indicate  that  the  airframe  shell 
structure,  which  accounts  for  approximately  28  percent  of  the  total  substruc- 
ture airframe  weight,  did  not  deform  appreciably,  due  in  part  to  the  stabiliz- 
ing effect  of  the  end  closures  (tension  rods),  during  the  27.5  ft/sec  vertical 
Impact.  Consequently,  the  occupant  responses  may  not  be  altered  very  much  from 
that  which  would  be  expected  if  the  upper  shell  structure  flexibility  were 
ignored  and  the  mass  and  inertia  effects  were  accounted  for.  This  situation  is 
particularly  significant  since  the  32  mass,  57  member  model  and  test  results 
show  that  the  shell  structure  motion  does  not  pose  any  lethal  threat  to  the 
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Figure  4-6.  Floor  and  occupant  peak  responses  as  a 
function  of  analysis  filter  cutoff 
frequency. 


occupants  either  through  failure  or  excessive  motion.  If  prior  to  or  during 
the  analysis,  the  user  can  establish  that  the  representation  of  the  shell 
structure  is  not  critical,  insofar  as  the  occupant  and  floor  resnonses  are 
concerned,  a smaller  more  economical  math  model  can  be  pursued.  To  assess  the 
tradeoff  between  accuracy  and  model  analysis  cost  for  the  substructure  and 
impact  condition  described  in  this  report  the  following  KRASH  math  models 
were  investigated: 

• 16  mass.  32  member  symmetrical  half  structure  (runmod-1) 

• 6 mass,  fl  member  symmetrical  half  structure  (runmod-1) 

• 5 mass,  3 member  full  structure  (runmod-O) 

The  16  mass,  32  member  model  is  shown  in  Figure  4-7.  The  model  is  the 
same  as  the  base  32  mass,  57  member  model  except  that  the  mass  associated  with 
the  shell  structure  is  distributed  among  the  floor  masses.  The  6 mass. 
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8 member  model,  shown  In  Figure  4-8,  lumps  half  the  airframe  weight  at  mass 
locations  1 and  2.  The  stiffness  of  the  seat  Is  represented  by  2 members 
Instead  of  4 as  in  the  larger  16  mass  and  32  mass  models.  The  seat  cushion, 
seat  pan,  occupant  and  DR1  are  represented  In  the  same  manner  as  in  the  larger 
models.  The  8 mass,  16  mass  and  32  mass  model  analyses  are  performed  as  sym- 
metrical half-structure  symmetrical  Impact  conditions.  The  5 mass,  5 member 
full  math  model,  shown  in  Figure  4-9,  treats  the  substructure  and  occupant  as 
a series  of  axial  members. 

The  results  of  the  model  size  variation  investigation  are  presented  in 
Table  4-1,  along  with  the  test  results.  The  average  of  the  floor  peak  re- 
sponses obtained  from  the  analyses  are  shown  to  range  within  0.5  to  13.7  per- 
cent of  the  average  of  the  test  peak  values.  The  analytical  occupant  lower 
torso  responses  are  within  19  to  22  percent  of  the  test  value  for  the  first 
peak  and  between  3 and  7 percent  of  the  test  values  for  the  second  peak.  The 
primary  upper  torso  peak  response  obtained  by  analyses  are  from  15  to  24  per- 
cent higher  than  the  corresponding  test  value.  As  is  anticipated  the  smallest 
model  is  nearly  1/8  as  expensive  to  run  as  the  largest  model.  The  6 and  16 
mass  models  are  approximately  20  and  75  percent,  respectively,  as  expensive  to 
run  as  the  32  mass  model.  The  results  indicate  that  for  this  particular  situa- 
tion the  use  of  smaller  KRASH  math  models  can  be  used  to  assess  trends  with 
sufficient  accuracy  and  substantially  lower  costs. 
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Figure  4-8,  6 mass,  8 member  symmetrical  math  model, 
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Figure  4- 9,  S mass,  5 member  full  math  model. 


LOCKHEID 


4-14 


TABLE  4-1.  COMPARISON  OF  RESULTS  FOR  DIFFERENT  SIZE  MATH  MODELS (a) 
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SECTION  5 
CONCLUSIONS 


The  analysis  of  a twin-engine,  low-wing  airplane  substructure  subjected 
to  a 27.5  ft /sec  vertical  impact  has  demonstrated  program  KRASH' s capability 

to  quantitatively  represent  the  significant  dynamic  response  phenomena,  namely: 

• Primary  floor  acceleration  magnitudes  and  times  of  occurrence 

• Occupant  response  magnitudes,  time  of  occurrence  and  response  shape, 
particularly  the  "camel-hump"  effect 

• Occupant  and  structure  motion  trends 

The  results  of  the  sensitivity  investigation  using  program  KRASH  indi- 
cate that  for  the  parameters  varied,  the  structure,  and  the  impact  condition 
evaluated,  user-selected  data  for  input  into  KRASH  can  vary  as  much  as  20  per- 
cent and  still  provide  a reasonable  assessment  of  overall  dynamic  behavior. 
Furthermore,  the  range  of  dynamic  response  peak  values  obtained  from  variations 
in  KRASH  user  selected  input  data  is  comparable  to  the  spread  in  the  measured 
test  data  between  the  left  and  right  sides. 

The  results  of  the  test  and  analysis  correlation  and  sensitivity  studies 
provide  valuable  information  which  can  be  used  to  enhance  future  modeling  of 
crash  impact  conditions. 

For  some  structural  configurations  and  Impact  conditions  simple  approxi- 
mate models  are  a cost-effective  method  of  representing  large  structural  seg- 
ments with  acceptable  accuracy  for  qualitatively  assessing  dynamic  behavior 
and  response  trends. 
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MASS  REPRESENTATIONS  MEMBER  REPRESENTATIONS 


MASS  MOMENT  OF  INERTIA  OF  50-PERCENTILE  MAN,  161.5  lb  (Reference  6) 
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SEAT  BELT 

Seat  Belt  Attaches  from  Mass  30  to  Mass  11  and  from  Mass  30  to  Mass 


Estimated  seat  cushion  stiffness  “ 500  Ib/in 
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EXTERNAL  SPRING  REPRESENTATION 


MASSES  1-6 
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MASSES  7-12 


MASSES  14,  17,  18 


MASSES  20.  21  MASSES  24-28 


MASS  DISTRIBUTION  BASED  ON  VOLUME 


TABLE  V-1.  - MASS  LOCATIONS  AND  PROPERTIES 
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